Overlap of quasiparticle random-phase approximation states for nuclear matrix 
elements of the neutrino-less double beta decay 
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Quasiparticle random-phase approximation (QRPA) is applied to two nuclei, and overlap of the 
QRPA excited states based on the different nuclei is calculated. The aim is to calculate the overlap 
of intermediate nuclear states of the double-beta decay. We use the like-particle QRPA after the 
closure approximation is applied to the nuclear matrix elements. The overlap is calculated rigorously 
by making use of the explicit equation of the QRPA ground state. The formulation of the overlap 
is shown, and a test calculation is performed. The effectiveness of the truncations used is shown. 
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Needless to say, neutrino physics is important for par- 
ticle physics as it provides information on lepton-number 
violation, Majorana nature of neutrino, and neutrino 
mass, see e.g. P, Q, which are the aspects beyond the 
Standard Model. Neutrino physics is also very interest- 
ing for nuclear physics |3l-[l0l|. because an accurate cal- 
culation of nuclear matrix elements of the neutrino-less 
double-beta {0v(3f3) decay is necessary, e.g. 0, [l^l, along 
with the experimental half-life value of the initial nu- 
cleus in order to determine the neutrino mass in one of 
few methods. This calculation is a challenging opportu- 
nity to test the capability of the techniques of theoretical 
nuclear physics for calculating many-body correlations in 
heavy nuclei. 

One of the important theoretical methods to calcu- 
late the nuclear matrix elements is the proton-neutron 
quasiparticle random-phase approximation (pn-QRPA), 
a method that has been improved significantly in the past 
few decades. A few of the milestones include the cor- 
rect treatment of the effects of the Pauli-exclusion prin- 
ciple in the intermediate states by the renormalizcd pn- 
QRPA [llj, fulfilment of the Ikeda sum rule by the fully 
renormalizcd pn-QRPA [12[, inclusion of the pn-pairing 
correlations [1^, and extension to the deformed states 
[l^ . Nevertheless, the nuclar matrix elements calculated 
by various approaches including the pn- QRPA arc dis- 
tributed in a range of a factor of 2 (3l.1T5|. 

Wc calculate the nuclear matrix elements of the 0^/30 
decay by employing the likc-particle QRPA [igI . [17| . 
which can be applied after the closure approximation 
0, [IS Il8l420j is used. The advantage of this approach 
is that the intermediate states are free from the Pauli- 
exclusion-principle problem without any modification be- 
cause they are states of even-even nuclei. In addition, we 
calculate the overlap of the intermediate states obtained 
by the like-particle QRPA wit h g reater accuracy than 
that in the previous studies [3, [3 [2l| - |23| . It is empha- 
sized that the importance of the overlap of the interme- 
diate states has been pointed out in Ref. [l3| in terms of 
deformation. The purpose of this study is to demonstrate 
the feasibility of the overlap calculation using the ground- 
state wave function of the like-particle QRPA explicitly. 
The equation of the QRPA ground state has been known 



for decades, e.g. [24|. To the best of our knowledge, how- 
ever, this is the first time that a numerical calculation 
has been carried out rigorously. 

The axial and parity symmetries of the nuclei are as- 
sumed throughout this paper. The z-component of the 
angular momentum is denoted by for nucleon state a 
and by Km for nuclear state m. tTq, and tt^ are used to 
denote the parity. Hereafter, we refer to the like-particle 
QRPA as simply QRPA. After the closure approximation 
is applied to the nuclear matrix elements of the de- 
cay, e.g. Refs. [3, [13, [Hi, one of those matrix elements, 
the Gamow- Teller type as an example, is written as 

M^T = {F\Y,h+{n,.Ea)n{i)-cT{j)T+{i)T+{j)\I) 

= 5IE<"">+(^12,^a)T(l)-^(2) 
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where \F) and |/) denote the final and initial nuclear 
states of the decay, and the ground states of the QRPA 
are used. h^{rij,Ea) is the neutrino potential [2| with 
fij = ki ~ '"jl, and Ea being the average energy of the 
intermediate nuclear states, i (j) indicates a nucleon. 
(T denotes the spin-Pauli matrix, r-|_ is the raising op- 
erator of the z-component of the isospin. An arbitrary 
single-particle basis {a} is introduced, and the creation 
and annihilation operators of the single-particle state are 
denoted by and Cq , respectively. The creation and an- 
nihilation operators of the excited state m of the QRPA 
based on the initial state are denoted by and O^, 
respectively, and those based on the final state are de- 
noted by and O^, respectively. The completeness 
using these operators is used in Eq. ([1]). There is no se- 
lection rule for the intermediate states with respect to 
(KmT^m)- Simplified notations K = Km and tt = tt™ are 
used hereafter. 

We express |/) in the form [23| 
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where \i) is the HFB ground state of the same nucleus 
as the one |/) describes, and Vj^ is a generator of the 
QRPA ground state. Afi is the normahzation factor. In 
this paper, for any equation only referring to the initial 
state denoted by / or i, we have also provided the corre- 
sponding equation referring to F or f. Latter equations 
are omitted. We have [0(,l,Of^,] = in the QRPA, if 

(Kn) 7^ {Km'T^m')', hence, v^/^ ''s with different (K'tt') 
are separately determined by 



0. 



(3) 



A general quasiparticle basis {^} based on the initial 
state is introduced by a^|i) = 0, where the is an- 
nihilation operator, /i ~ (9^, tt^, J^, i^) being the label 
of a general quasiparticle state, denotes a proton or 
neutron, and is the label specifying a state in the sub- 
space (g^, TT^, j^). Notation — /i is used for expressing 
(gu. TTu, —7m, iu.)- We use the canonical-quasiparticle ba- 
sis [131 {m} for efficiency of the QRPA calculation [26| . 

'A 



V} 'is expressed as 



ajj^al^ and a^Ja^J in Ec^. are the fermion image of 

the boson [U; a condition is introduced that C^^JlJ/ 
does not vanish, only if + = K' , f^, + j^, = -K', 
and TT^TTjy = TT^'TTiy/ = tt'. Wc order the canonical- 
quasiparticle states and restrict ^ < ly, fi' < ly' in 
n'l' without loss of generality. 



The solution of the QRPA equation gives us 



where -I- = Km' and tt^tt^ 
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We define matrices 



, (6) 



(7) 



where the QRPA solutions having {K'tt') are used. The 
negative integers of the index correspond to Matri- 
ces C'-K'''')P, x'-^''''^'^, and F('f<"'^')^ are also 
introduced in the same way. 

(j{K'^')i ig obtained, by 
ignoring the exchange terms (the quasi-boson approxi- 
mation), as follows: 



1 
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where suffix T stands for transpose. Practically, 

a singularity. 



The relation between the two HFB states can be writ- 
ten as 



N) = ;^exp 
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/ is the unit matrix of the size of matrix I?, which is 
defined by 



(11) 



riT is the dimension of the subspace (g^, tt^, j^). Z^^j^ 
is not equal to only for those ^ and v that satisfy 
jfi+jl = and TT^TTy = We restrict > in Eq. 
The unitary transformation from basis {a^^, a^^} to ba- 
sis {a^^ , o-^-^j} is given by 



E rr^IFl Ft I T,/F2 F 



(12) 



and its Hermite conjugate equation for — /i with = j^, 
and TT,, = tTm/. Tlf"} and T/^^, can be obtained from 
the volume integral of the product of the canonical- 
quasiparticle wave functions |27j of the two bases. D^-y 
is given by 



D = -( -^T'^^ 



(13) 



r^^).-T;l/,^^,J-l,---,nT. (14) 

Matrix T^^'^ is defined in the same way as matrix D. 
Practically, again, does not have a singularity. 

Now, we expand and truncate the overlap matrix ele- 
ment with respect to v^J^ ' and ' as 
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We test up to the second-order terms G^^,, which use 



both v'-p^''^' and f^'^^'^'-'^but only with (i^iTTi) = [K2t:2) 



vf''''^ and f^^^"^) 
in Eq. (fTS)) (actually G^{^, is negligible in most of the 
overlap matrix elements, as shown later). Up to the 
fourth-order terms are included in normalization factors 
Mi and Nf, because its convergence of the {)-expansion 
is slow as compared to the un-normalized overlap ma- 
trix elements (which are the result of the numerical 
test). Equations (fT6|) — (fTS]) can be calculated using Xjj^ , 
^'^'"'^^ those of F, T.^^i ^^^2 



and Dfj^iy. The 



concrete equations will be given in the forthcoming full 
paper. 

We use the code of the HFB approximation and 
that of the QRPA developed by us [1^ . The wave func- 
tions are treated, in both the codes, on a mesh within 
the cylindrical box and are discretized by the vanish- 
ing boundary condition at the edge of the box. The 
HFB equation is solved using a cutoff of the quasipar- 
ticle energy at 20 McV for convenience in performing the 
tests. We transform the wave functions of the quasipar- 
ticle states to those of the canonical-quasiparticle states 
before the QRPA calculation 

2^Mg and ^^Si are used in this paper for \i) (|/)) and |/) 
(|F)), respectively, with the Skyrme parameter set SkM* 
|30| and the volume pairing density functional [3l|. The 
properties of the HFB ground states are shown in TableU 
The total dimension of the HFB space is ~330, including 
those with negative jz- Strength of the volume pairing 
Gn = —270 MeV fm'^ is used for the neutrons of both 
26Mg and 26Si, and Gp = -150.0 and -270.0 MeV fm^ 
is used for the protons of ^^Mg and ^^Si, respectively. 
These values were chosen so that the mean fields of the 
ground states arc similar between the two nuclei. 



TABLE I: Properties of the HFB ground states of ^''Mg and 
^®Si used in this paper. /3p and Ap denote, respectively, the 
quadrupole deformation and the averaged pairing gap of the 
protons. /3n and An denote the same for the neutrons. 



nucleus 


/3p 


Ap (MeV) 


/3n 


An (MeV) 


26Mg 


-0.199 


0.794 


-0.195 


1.510 


26Si 


-0.224 


0.865 


-0.206 


1.402 



We show the results of the calculation of [Kn) = (0+) 
below. Let 'Tip and *Jl/ be the number of the two- 
canonical-quasiparticle states associated with \F) and 
|/), truncated by the cutoff occupation probability for 
calculating Eqs. (fT7|) and (fT8|) (those with larger occu- 
pation probabilities than the cutoff arc used). This is 
another truncation after the 20-McV cutoff. The con- 
vergence of the overlap matrix elements is shown with 
respect to OTf in Fig.[TJ The same value of the cut- 
off is applied for the two bases, and we have ^Rf — 91/. 
It is seen that '^Ip + = 350 is sufficient for the con- 
vergence. The total number without the truncation is 
~3300; thus, this truncation is rather efficient. |/) and 
\F) have different configurations at the Fermi surface; 
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FIG. 1: (Color online) Ten diagonal overlap matrix elements 
having the largest absolute values, as functions of Vip + 01/. 
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FIG. 2: (Color online) Twenty largest absolute values of 
the diagonal overlap matrix elements. Those up to the sec- 



ond order with respect to t)}"^' 



and i) 



(0+) 



are shown in the 



descending order. (T^itti) 7^ (0+) are not included in 
([T7|) and G^^, 91f + 01/ = 134 was used; see Fig.[Il 



therefore, the high-energy components of 0,^J, and 
leaving the configuration around the Fermi surface intact 
do not contribute to the overlap matrix elements. On the 
other hand, it is necessary to calculate A// ([T9|) and Mp 
without this truncation. 

The major diagonal overlap matrix elements are shown 
in Fig. [21 obtained with 91f + 51/ = 134. It is observed 
that the contribution of G^^^, ^8]) is negligible and that 

^mm' C3) is not significant for the small matrix el- 
ements. G^^i (|16p is sufficient in most of the matrix 
elements omitted in that figure. 

The contribution of (Xitti) 7^ (0+) to the major over- 
lap matrix elements through G^^^ PT|) is shown in Fig. [31 
calculated with OT/r -|- 01/ = 350 and max \K\ = 4. We 
also calculated the contribution of {KiHi) ~ (0—) and 
(1—) and found that it was smaller than that of the 
positive parity by at least an order of magnitude; thus. 
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FIG. 3: (Color online) Twenty largest absolute values of 
the diagonal matrix elements of the overlap. G^^/ (|17|l was 
calculated with (Kini) = (0+) — (4+) (all terms), and the 
terms with only (0+) {Ktv = 0+ term). ^Xf + Vli = 350 is 
used. G^^; is not included. QRPA solution number m does 
not necessarily correspond to that of Fig. [2] 

only the positive parity is used. The contribution of 
(KiTTi) 7^ (0+) is very small to all of G^^^ except for 
TO = 1, which is one of the spurious states. Actually, our 
method should be applied only to the cases that do not 
have large fluctuations of the particle number so that the 
spurious states are not crucial to the nuclear matrix ele- 
ments (IlJ. Mp and Afi require \K\ of up to 3 with both 
parities. 

In summary, the overlap matrix elements of the QRPA 



states based on the ground states of different nuclei have 
been calculated using the QRPA ground states explicitly 
for relatively light nuclei with the Skyrme and the con- 
tact volume pairing energy functionals. The most im- 
portant finding of this study is that the bold truncations 
are allowed in the calculation of the un-normalized over- 
lap matrix elements. The normalization factors need a 
less-truncated calculation; however, the amount of this 
calculation is reduced tremendously by identifying |/) 
and \i) in each factor. Considering this advantage and 
the performance of the modern parallel computers, we 
believe that there is no reason to avoid the explicit wave 
function of the QRPA ground states. 
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